CORRELATED CONTINUOUS TIME RANDOM WALKS 
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Abstract. Continuous time random walks impose a random waiting time before 
each particle jump. Scaling limits of heavy tailed continuous time random walks 
are governed by fractional evolution equations. Space-fractional derivatives describe 
heavy tailed jumps, and the time-fractional version codes heavy tailed waiting times. 
This paper develops scaling limits and governing equations in the case of correlated 
jumps. For long-range dependent jumps, this leads to fractional Brownian motion 
or linear fractional stable motion, with the time parameter replaced by an inverse 
stable subordinator in the case of heavy tailed waiting times. These scaling limits 
provide an interesting class of non-Markovian, non-Gaussian self-similar processes. 



1. Introduction 

Continuous time random walks (CTRW) separate IID particle jumps {¥„) by IID 
waiting times (Jn). CTRW models are important in applications to geology [7j, 
physics [28] , and finance [32] . In the case of heavy tailed waiting times, CTRW scaling 
limits are subordinated processes that are self-similar but non-Markovian [27]. Their 
transition densities are governed by fractional diffusion equations [25]. Fractional 
diffusion equations replace the usual integer order derivatives in the diffusion equation 
by their fractional-order analogues [29l[30]. Just as the diffusion equation dtu = ad^u 
governs the scaling limit of a simple random walk, the fractional diffusion equation 
d^u = ad°u governs the scaling limit of a CTRW with heavy tail jumps P(F„ > r) ~ 
r"" for < a < 2 and waiting times P(Jn > t) ~ for < /5 < 1. 

This paper develops limit theorems and governing equations for CTRW with cor- 
related jumps Yn = CjZn-j, where are IID and (c„) are real numbers (see 
Section 2 for precise conditions). These CTRW models are useful for correlated 
observations separated by random waiting times, which are common, for example, 
in finance [31]. Scaling limits of the partial sum process S(t) = Yi + ■ ■ ■ + Y^t] in 
the case of long range dependence include fractional Brownian motion (FBM) for 
light-tailed jumps [TBI [55] , and linear fractional stable motion (LFSM) for heavy- 
tailed jumps [H [211 ES]- Letting Tn = Ji + ■ ■ ■ + Jn the time of the nth jump, and 
Nf = max{n : T„ < t} the number of jumps by time t > 0, the scaling limit of the 
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CTRW S{Nt) is a FBM or LFSM subordinated to an inverse stable subordinator, 
which is connected to the local time of a strictly stable Levy process [2Bj, or the 
supremum process of a spectrally negative stable Levy process [11]. This extends 
the results of [271 E] to the case of dependent jumps. We also discuss some interest- 
ing properties of these self-similar limit processes, and governing equations for their 
probability densities. 



2. Results 

Let {Zn, —oo < n < oo} denote IID random variables that belong to the strict 
domain of attraction of some stable law A with index < a < 2. This means that 
the sequence of partial sums P{n) = Zi + ■ ■ ■ + Zn satisfies anP{n) =^ A for some 
an > 0, see Feller [HI p.312-313] or Whitt [36l p.114-115]. 

The particle jumps that we consider in this paper are given by the stationary linear 
process — oo < n < oo} defined by F„ = Yl'jLo '^j^n-j, where Cj are real constants 
such that l^jl^ ^ ^ some p G (0,a). This condition ensures that the the 

series Yl'jLo'^j^n-j converges in Lp{F) and almost surely (see Avram and Taqqu [2]). 

Let J„ > be IID waiting times, T„ = Ji + - ■ ■ + J„ the time of the nth particle jump, 
and Nt = max{n : T„ < t} the number of jumps by time t > 0. Let S{n) = Yi+- ■ ■+Yn 
denote the location of the particle after n jumps, so that the continuous time random 
walk (CTRW) S{Nt) gives the location of the particle at time t > 0. Suppose that 
Jn belongs to the domain of attraction of some stable law D with index < (3 < 1 
and D > almost surely. Hence 6„T„ =^ D for some norming constants 6„ > 0. Let 
b{t) = and take b{t) an asymptotic inverse of the regularly varying function l/b{t), 
so that tb{b{t)) 1 as t ^ oo [27J. 

Let {A{t),t > 0} and {D{t),t > 0} be stable Levy processes with A{1) = A, D{1) = 
D, respectively. Note that {D{t),t > 0} is a stable subordinator of index j3, hence its 
sample functions are almost surely strictly increasing [8^ p. 75]. Therefore, the inverse 
or hitting time process of {D(t),t > 0}, 

Et = M{x > : D{x) > t}, Vt > 

is well defined and the function t ^ Et is strictly increasing almost surely. 

Our first result shows that the CTRW scaling limit in the case of short-range 
dependence is quite similar to the case of independent jumps studied by Meerschaert 
and SchefHer [27j. 

Theorem 2.1. Under the conditions of this section, suppose that < a < 2, Cj > 
and Cj<oo for some p G (0, a) with p < 1, and that one of the following holds: 

(a) < a < 1; or 

(b) Cj = for all but finitely many j ; or 

(c) 1 < a < 2, Cj is monotone and < oo for some p < 1. 
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Then we have 

(2.1) w~\^,^^SiN,t) ^ AiE,) 

05 c — > oo in the Mi topology on D{[0, oo),M), where w = J2j ^j- 

In view of Avram and Taqqu [2j , the convergence in (12. ip cannot be strengthened to 
the Ji topology. Note that the processes {A{t)^t > 0} and {Et, t > 0} are self-similar, 
that is, for every constant c > 

{A{ct),t> 0} = {ci/"A(t),t > 0} 

and 

{^rf,t>0} = {c^Ei,t>0}, 

where = means equality in all finite dimensional distributions. It follows immediately 
that the scaling limit {A{Et),t > 0} in (12.11) is self-similar with index f3/a. When 
< /3 < 1/2, the inner process Et in (12. ip is also the local time at zero of a strictly 
stable Levy motion [26]. When 1/2 < f3 < 1, the inner process Et is also the 
supremum process of a stable Levy motion with index 1/(3 and no negative jumps 

m 

Let d^git) denote the Caputo fractional derivative, the inverse Laplace transform 
of s^g{s) — s^~^g{Q) where g{s) = e~** g(t) dt is the usual Laplace transform of g. 
Let d^^f{x) denote the Liouville fractional derivative, the inverse Fourier transform 
of {±ik)°'f{k), where f{k) = f^^e~'^^^f{x)ds is the usual Fourier transform. The 
stable random variable A{t) has a smooth density with Fourier transform e~^'^'^^'^ 
where ip{k) = a[p{ik)°' + q{—ik)°'] with < p,q < I and p + q = 1 [23]. Then the 
limit A{Et) in (12. ip has a density h{x, t) that solves the fractional diffusion equation 
d^h = apd^h + aqd^^h, see |27j. 

Next we consider the CTRW scaling limit for heavy-tailed particle jumps with 
long-range dependence. To simplify the presentation, we assume = n~^/" (domain 
of normal attraction) and power-law weights; namely Cj ~ Coj^~^~^^°' as j — >■ cxo, 
for some Cq > 0. The case < H < 1/a means the stationary sequence {Yn} 
has short-range dependence, while the case 1/a < H < 1 means {Yn} has long- 
range dependence. The scaling limit of CTRW with short-range dependence has 
been partially covered by Theorem 12. 1[ The rest of the cases are treated in Theorems 
[221 and O below. 

We will make use of the following definition. Given constants a G (0,2) and 
H G (0, 1), the a-stable process {Lo,jj{t),t G M} defined by 

(2.2) L^,H{t) = [ \{t- s);~'/" - (-5);-'/"] A{ds) 

is called a linear fractional stable motion (LFSM) with indices a and H. In the 
above, a+ = max{0,a} for all a G M and {A{t),t G M} is a two-sided strictly stable 
Levy process of index a with A{1) = A given at the beginning of Section 2. Namely, 
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n a.s n oo. Because of this, {La,H{t),t G M} defined by (2.2) differs 

from the LFSM in Theorem 4.7.2 in [36j by a constant factor. Note that, when 
H = 1/a, La^nit) = A{t) for all t > 0. When H ^ '^/a-, the stochastic integral in 
fl2.2l) is well-defined because 



dr < CO. 



See [33j, Chapter 3]. 

By (12.21) . it can be verified that {La^H{t),t G M} is if-self-similar with stationary 
increments \S3\ Proposition 7.4.2]. It is an a-stable analogue of fractional Brownian 
motion and its probabilistic and statistical properties have been investigated by many 
authors. In particular, it is known that 

(i) If 1/a < if < 1 (this is possible only when 1 < a < 2), then the sample 
function of {La^H(t),t G M} is almost surely continuous. 

(ii) If < H < 1/a, then the the sample function of {La,H(t),t G M} is almost 
surely unbounded on every interval of positive length. 

We refer to [33, Chapters 10 and 12] for more information. 

Theorem 2.2. We assume the setting of this section. If 1 < a < 2, 1/a < H < 1, 

and Cj ~ Coj^^^^^^" as j — * oo for some Cq > 0, then as c ^ oo 

(2.3) [b{c)]-''S{N,t) L^,H{Et) 

in the Ji topology on D([0, oo), M), where Ki = Coa/{Ha — 1). 

The topology on D([0, oo),M) in Theorem 12.21 is stronger than that in Theorem 
12.11 thanks to the fact that La^nit) is a.s. continuous whenever 1/a < H < 1. 

Observe that the case when < H < 1/a and the constants Cj {j > 0) are not 
all nonnegative is left uncovered by Theorems 12.11 and 12.21 Because of Property 
(ii) of {La,H(t),t G M}, the limiting process does not belong to the function space 
Z)([0, oo), M). Nevertheless, we have the following theorem. 

Theorem 2.3. We assume the setting of this section. If < a < 2, < H < 1/a, 
Cj ~ Coj^^^^^^"' as j — > oo for some cq > 0, and J2j Cj = 0, then 

(2.4) [6(c)]-^5(iV,,) ^ L^MEt) 

as c ^ oo, where means convergence of all finite- dimensional distributions and 
Ki = coa/{Ha — 1). 

It is interesting to note that the constants in Theorems 12.21 and 12.31 are determined 
by Co, a and H in the same way. But Ki is positive is when 1/a < H < 1, and is 
negative when < H < 1/a. 

It follows from the self-similarity of {Lc,,H{t),t G M} and {E{t),t > 0} that the 
scaling limits in (12. 3p and (12. 4p are self-similar with index H(3. When 1/a < H < 1, 
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it can be seen that {La^H{Et),t > 0} has continuous sample functions almost surely. 
However, if < H < 1/a, then {La,H{Et),t > 0} is almost surely unbounded on every 
interval of positive length. It would be interesting to further study the properties of 
the process {La,H{Et),t > 0}. 

We mention that both Theorems 12.21 and 12.31 can be extended to (Z„) in the strict 
domain of attraction of A and (cj) regularly varying at oo with index H — 1 — 1/a, 
using a slightly different normalization in (12.31) depending on and the probability 
tail of Zn, compare [IJ. 

Finally we consider the case a = 2. If {A{t),t G M} in fl2.2p is replaced by 
ordinary two-sided Brownian motion, then fl2.2l) defines a fractional Brownian motion 
Wh = {WH{t),t G M} on M of index H, which is a Gaussian process with mean zero 
and covariance function 

E[WHit)WHis)] = ^ 0^1''' + i^i'"" - 1^ - ■ 

Theorem 12.41 gives the CTRW scaling limit for light-tailed particle jumps with 
long-range dependence. 

Theorem 2.4. We assume the setting of this section. If a = 2, E[Z„] = 0, IE[^^] < 
J2j ^'j < Var[S'(n)] = varies regularly at oo with index 2H for some 
< H <1, andE[S{nfP] < K2 \E{S{nf)Y for some constants K 2 > and p > 1/H, 
then as c 00 

(2.5) %U]SiN,t) ^ WniEt) 

in the Ji topology on D([0, 00), M). 

Note that it is not difficult to provide examples of sequences of IID random variables 
{Zn} and real numbers {cj} that satisfy the conditions of Theorem 12. 4^ see [121111]. It 
follows from the results of Taqqu [35] that the conclusion of Theorem 12.41 still holds if 
the linear process {Yn} is replaced by the stationary sequence {g{^n)}, where {^„} is 
a stationary Gaussian sequence with mean 0, variance 1 and long-range dependence, 
and g G L^(e~^' ^"^dx) is a function with Hermite rank 1. 

Theorem 12 . 41 contains the case H = 1/2 where Wnit) = A{t) is a standard Brownian 
motion. This includes the situation of mean zero finite variance particle jumps, and 
heavy tailed waiting times between jumps. In this case, the CTRW scaling limit A{Et) 
has a density h{x, t) that solves the time-fractional diffusion equation d^h = a d^h, 

see [27]. Since {WH{ct),t > 0} = {c^WH{t),t > 0}, the scaling limit in ([23]) is self- 
similar with index H(3. Some results on large deviation and sample path regularity 
have recently been obtained for {WH{Et),t > 0} in [26] . 

In the case of finite mean waiting times, the CTRW scaling limit is essentially the 
same as for the underlying random walk. If /i = EJ^ < 00, then fiNt/t — >■ 1 almost 
surely as t 00, and a simple argument along the lines of the proof of Theorem 12.11 
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shows that w^^a[c]S{Nct) =^ A{t/fj,) in the Mi topology on /^([O, oo), M). Theorems 
12.21 12.31 and 12.41 can be extended similarly. 

An easy argument with Fourier transforms shows that the density h{x, t) of La^ui^) 
solves dth = aHt'^^~^[apd'^h + aqd'^^h]. A similar argument shows that the density 
of Wnit) solves dth = 2Ht^^~^a&lh. An interesting open question is to establish 
the governing equation for the CTRW scaling limit in (12.31) and (12.51) . This is not as 
simple as replacing the first time derivative by a fractional derivative in the governing 
equation for the outer process, since the t variable also appears on the right-hand 
side, so that Theorem 3.1 of [3] does not apply. 

3. Proofs 

The proofs in this section are based on invariance principles for stationary sequences 
with short or long-range dependence (see, for example, [36]) and the CTRW limit 
theory developed in [27j. Due to the non-Markovian nature of the CTRW scaling 
limits in this paper, standard subordination methods can not be applied directly. 
Instead we apply continuous mapping-type arguments to prove Theorems 12. ![ 12.21 
and 12.41 The proof of Theorem 12.31 is quite different and relies on a criterion for the 
convergence of all finite-dimensional distributions of composite processes established 
by Becker-Kern, Meerschaert and SchefHer |6]. 

Recall that J„ > are IID waiting times, = Ji + ■ ■ ■ + Jn the time of the nth 
particle jump, and Nt = max{?T, : T„ < t} the number of jumps by time t > 0. Since 
Jn belongs to the domain of attraction of some stable law D with index < (3 < 1 and 
D > almost surely, with 6„T„ =^ D for some norming constants 6„ > 0, the sequence 
bn varies regularly with index — see [18]. Then the asymptotic inverse b{t) of 1/6 
varies regularly with index f3, see Seneta [31]. Recall that the stable Levy motion 
{D{x),x > 0} with -D(l) = D is a. stable subordinator of index /3 and thus is strictly 
increasing [SI p. 75]. Its inverse or hitting time process Et = inf{x > : D{x) > t} 
is strictly increasing with continuous sample paths, has moments of all orders, and 
its increments are neither stationary nor independent [27]. Bingham [10] shows that 
Ef has a Mittag-LefHer distribution, and gives a differential equation that governs its 
finite dimensional distributions. 

Proof of Theorem \2.1\ Corollary 3.4 in |27] shows that b{c)~^Nct =^ Et as c —>■ oo 
in the Ji topology on D{[0,oo), [0, cx))). Note that 6(c) — > oo as c ^ oo since this 
function is regularly varying at oo with index f3 > 0. Theorem 4.7.1 in Whitt 
shows that anS{nt) =^ wA(t) in the Mi topology on D([0, oo),M). Since the Ji 
topology is stronger, and since the waiting times (J„) are independent of (Y„), we 
have 

(«[6(c)]^(Kc)t),6(c)-iAr,,) ^ iAit),Et) 

in the Mi topology of the product space -D([0, oo),]R x [0, oo)). Of course, this last 
statement also follows from Theorem 3.2 in [9j. 
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Since the process {E^, t > 0} is strictly increasing and continuous, Theorem 13.2.4 



Proof of Theorem \2.2. Recall that b{c)~^Nct =^ Et in the Ji topology on D([0,oo), 
[0, oo)) [271 Corollary 3.4]. Theorem 4.7.2 in Whitt [36], originally due to Astrauskas 
[T], shows that n~^S{nt) =^ KiLa^nit) in the Ji topology on D([0, cxd), M), where 
Ki = co(y/{Ha — 1). 

Since {Nt,t > 0} is independent of {S{n),n > 1}, we have 



in the product space. Combining this with [3B1 Theorem 13.2.4] yields fl2.5l) in the Mi 
topology. Since both processes {La,H(t),t > 0} and {Et,t > 0} are continuous, and 
the latter is strictly increasing, one can apply Theorem 13.3.1 in [36] to strengthen 
to the conclusion to convergence in the Ji topology. This proves Theorem 12. 2[ □ 

Proof of Theore'm \2.3\ . It is sufficient to show that for all integers m > 1, < ti < 
■ ■ ■ < tm, we have 



as c ^ oo. For this purpose, we will make use of Proposition 4.1 in [6], which 
provides a useful criterion for the convergence of all finite-dimensional distributions 
of composite processes, and Corollary 3.3 in [2T] which is concerned with convergence 
of finite dimensional distributions of weighted partial sums of IID random variables. 

We will adopt some notation from [B]. For t = (ti, . . . , t^) and c > 0, let Pc := 
be the distribution of b{c)~^ {Nd-^, ■ ■ ■ ,Nct^), and let p := p* be the distribution of 
{Et„...,EtJ. Since b{c)''N,t Et in the Ji topology on D([0,oo), [0,oo)) [271 
Corollary 3.4], we have pc ^ P as c oo. 

It follows from the definition of {Yn} that, for every x >0, S{nx) can be rewritten 

as 



where Ck = Ck ii k > and Cfc = if A; < 0. Under the assumptions of Theorem 12.31 
we have E^-oo < oo, E^-oo^fc = and 




□ 



{[bic)]-''S{b{c)t),bic)-'N,t) (Ki L^,H{t),Et) 



(3.1) 



6(c)-^(5(iV,,J,...,^(iV,,„)) ^ Ki{L^^H{Et,),...,L^MEtJ) 



(3.2) 




oo 



oo 




k=n 



k=n 
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a.s n —>■ oo. Thus, the conditions of Theorem 5.2 in [21J are satisfied with ipln) = 
^H-i/a^ a = -Ki (recall that Ki = coa/{Ha - 1)), 6 = and A = 0. It follows that 

(3.3) S{nt) KiLa^nit) as n ^ oo. 

For any x = (xi, . . . , Xm) G M^, let /ic(x) be the distribution of b{c)~^ (^S {b{c)xi) , 
. . . , 5'(6(c)xm)) and let i^(x) be the distribution of Ki (Lq, jj(xi), . . . , Lq, ^^(xm))- Then 
for every c > 0, the mapping x i— /ic(x) is weakly measurable. Since the linear 
fractional stable motion {La^H(t),t > 0} is stochastically continuous, the mapping 
X I— s> z/(x) is weakly continuous. Moreover, it follows from (13.31) that, for every x G M™, 
/ic(x) =^ z/(x) as c — > cxD. 

As in [6j we apply a conditioning argument and the independence between the se- 
quences {Yn} and {Jn} to derive that the distribution of b{c)~^ [S{Ncti), ■ ■ ■ , S{Nct^)^ 
can be written as J^rn /^c(x) dpd'x.), which is a probability measure on M"*. Similarly, 

the distribution of the random vector Ki(^La,H{Et-^), ■ ■ ■ , La,H{Et^)^ can be written 
as J^m z^(x) (ip(x). 

Therefore, (13. ip follows from Proposition 4.1 in P| once we verify that, for every 
X G M™, /ic(x'''^^) =^ z/(x) for every sequence {x'^'^^} C M.^ that satisfies x'^'^^ ^ x as 

C — > CXD. 

The last statement is equivalent to 

(3.4) c-^ {S{cx['^), . . . , S{cxl^^)) {L^,h{xi), . . . , W(x^)) 

whenever x^'^^ ^ x as c ^ oo. This is stronger than (13. 3p . where the fixed time- 
instants < xi < X2 < ■ ■ • < Xm are replaced now by x^'^^ . . . , Xm . Our proof of (13. 4p 
is a modification of the proof of Theorem 5.2 in Kasahara and Maejima [21j . 
To this end, we define the step function r i-^ Ac{r) on M by 

Then for any function g on M, as in [211 p. 88], we define 

/oo , oo 

g{r)dA,ir) = -^ E 3& 

j = '00 

By using (13. 2p . (13. 5p and (13. 6p we can rewrite S{cx) (x > and c > 0) as 

c-^5(cx) = ^ Yl ^^r^i ^A^i 

j=—oo ^ k=l — j ' 



(3.7) 

1 
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where the integrand gc{x,r) is given by 

[cx] — [cr] 



(3.8) 

^fc=-[cr]+l fe=[cx]-[cr]+l 

In the above, we have used the fact that X]jl-oo^i = to derive the second equahty. 
It follows from (13.71) that (13.41) can be rewritten as 

(3.9) I / glxf,r)dAlr)\ ^ [ gix,,r)dAir) 



i=l JK J i=l 



where g{x, r) = {x — r)^^^^" — (— r)^^"*^^" is the function in (12. 2p . 

Now let us fix X = {xi, . . . , Xm) G and an arbitrary sequence {x*^^)} C 
that satisfies x*^"^) — > x as c ^ oo. By Corollary 3.3 in [21j (with /^(■) being taken as 
gc{xf^\ ■)), the convergence in (13. 9p will follow once we verify that for every 1 < i < m 
the following conditions are satisfied: 

(Al)' for c/r-almost every r G M, 

(3.10) 9c{xt\rc) — K,g{xi,r) 

whenever rc ^ r as c ^ oo. 
(A2)' for every T > 0, there exists a constant 3 > a such that 

.p/ i.(4V)r...)<=., 

c>l J\r\<T 

where Pc{r) = [cr]/c, and 
(A3)' there exists an e > such that 

(3.12) lim limsup / < |5'c(a;i'^\ ''^) + \9c{xf'\r)\'^~^'^ \ dpc{r) = 0. 

T— >oo c-^oo J\r\>T ^ ^ 

For simplicity of notation, we will from now on omit the subscript i. To verify 
Condition (Al)', note that by the property of {cfc}, we have 

1 ^ ~ \ -fsTir^-i/" if r > 0, 



1 _ f 

fc=[crl+l ^ 



if r < 0, 



and the convergence is uniform in r on every compact set in M\{0}. For any x G M+ 
and r G M, we may distinguish three cases r<0,0<r<x and r > x. By applying 
(I3.13P to (13. 8p we derive that, as c — oo, g^x, r) g{x, r) uniformly in (x, r) on every 
compact set in {(x, r) : x G M+, r G ffi\{0, x}}. This implies that gc{x^'^\rc) g{x, r) 
whenever r ^ x and — r as c — > oo. Hence (Al)' is satisfied. 
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To verify Condition (A2)', we take a constant {3 > a such that I3{H — 1/a) > —1 
and consider the integral 



(3.14) 



\<T 



1 



E 



fc=[cx] — [cr]+l 



1 



\j\<cT 



J2 

k=[cx]—j+l 



Let > 1 be a constant such that \ck\ < 2coA;^"^/"~^ for all k > N. We split the 
summation in the right-hand side of (13.141) according to whether [cx] — j < N or 
[cx] — j > N. Thanks to the fact that Cfc = we have 



(3.15) 



E 



1 



J2 

k=[cx] 



/3 



< 



,/3(H-l/a)+l 



O 



\j\<cT,[cx]-j<N 

for some constant K3 > 0. In the above we have also used the fact that there are at 
most + 1 non-zero terms in the summation in j. 
On the other hand, we have 



E 



\j\<cT,[cx]-j>N 



^P(H-l/a)+l 



(3.16) 



\j\<cT,[cx]-j>N 



J2 

k=[cx]—j+l 

{[x] - [j7c])'3(^-i/") 



/3 



\x - r|^(^-i/") dr 



\r\<T 



for some constants K4, K5 > and the last integral is convergent because (3{H — 
1/a) > —1. Combining fl3.15p and (13.161) yields 



(3.17) 



\r\<T 



< 



A'3 



E 



k=[cx] — [cr]+l 
J\r\<T 



dpcir) 

^_ ,/3(H-l/a) 



Since the first sum in (13.81) corresponds to x = 0, we see that (A2)' follows from 

(EIZD. 

The verification of (A3)' is similar to the above and hence is omitted. This finishes 
the proof of Theorem 12. 3[ □ 

Finally we prove Theorem 12. 4[ 

Proof of Theorem \2.4\ Recall that b{c)~^Nct =^ Et in the Ji topology [271 Corollary 
3.4]. Theorem 4.6.1 in Whitt p56j shows that, as n — 00, a~^S{nt) =^ Wnit) in the 
Ji topology on Z}([0, 00), M). This result is originally due to Davydov p6|, see also 
Giriatis et al. [HI p. 276]. Since the sequence (J„) is independent of (1^), we have 
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a -j S{b{c)t),b{c) ^Nct) =^ {Wnit), Et) in the product space, and then continuous 



mapping along with Theorem 13.3.1 in [36j yields (12.51) in the Ji topology. □ 



4. Discussion 



Self-similar processes arise naturally in limit theorems of random walks and other 
stochastic processes, and they have been applied to model various phenomena in a 
wide range of scientific areas including telecommunications, turbulence, image pro- 
cessing and finance [1^. The most prominent example is fractional Brownian motion 
(FBM). However, many real data sets are non-Gaussian, which motivates the develop- 
ment of alternative models. Many authors have constructed and investigated various 
classes of non-Gaussian self-similar processes. Samorodnitsky and Taqqu [33j pro- 
vide a systematic account on self-similar stable processes with stationary increments. 
Burdzy [T2l [T3] introduced iterated Brownian motion (IBM) which replaces the time 
parameter of a two-sided Brownian motion by an independent one- dimensional Brow- 
nian motion B = {Bt,t > 0}. In this paper we have shown that the limit processes 
of CTRWs with correlated jumps form a wide class of self-similar processes which are 
different from the existing ones. 

When < /3 < 1/2, the inner process Ef in (12.11) or (12.51) is also the local time at zero 
Lt of a stable Levy process, and the iterated process {WH{Lt),t > 0} is called a local 
time fractional Brownian motion (LTFBM) in [26], a self-similar process with index 
jSH and continuous sample paths. Large deviation and modulus of continuity results 
for LTFBM are developed in a companion paper |26]. Strassen-type law of the iterated 
logarithm has been proved by Csaki, Foldes and Revesz [15] for local time Brownian 
motion (LTBM, the case H = 1/2). It is interesting to note that our Theorem 12.41 
shows that "randomly-stopped stationary sequence" {(F„ : n < Nt),t > 0} belongs 
to the domain of attraction of {WH{Lt),t > 0} for all H G (0,1). This theorem 
provides a physical interpretation of the process {WH{Lt),t > 0}. 

One interesting property of LTBM is that its increments are uncorrelated (this 
follows by simple conditioning argument), but not independent. It has long been 
recognized that price returns are essentially uncorrelated, but not independent [5l [23] . 
Hence LTBM, the scaling limit of a CTRW with (weakly) correlated price jumps, may 
be useful to model financial price returns. This approach could provide an interesting 
alternative to the subordinated variance- Gamma model of Madan and Seneta [22l [T3] 
or the FATGBM model of Heyde [20]. 




LTBM has a close connection to fractional partial differential equations. Meer- 
schaert and SchefHer [27] and Baeumer and Meerschaert [3] showed that the proba- 
bility density u{x,t) of LTBM solves the fractional Cauchy problem. 




(4.1) 



dtu{t,x) = dlu{t,x). 
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Baeumer, Meerschaert and Nane [1] further showed that the density of the iterated 
Brownian motion solves the same equation (14. II) . As we mentioned at the end of 
Section 2, the connection between the hmit processes in this paper and fractional 
partial differential equations remains to be investigated. 
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